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EXECUTIVE  SUMMARY 


This  report  suggests  that  some  current  Navy  tasks  may  benefit  somewhat  if 
some  of  the  older  literature  is  re-examined  with  respect  to  the  current 
objectives  and  goals  of  these  ongoing  tasks.  This  report  uses  early  60 's 
mustard  .gas  Navy  paint  data  as  an  example  to  demonstrate  how  older  data  may 
still  remain  useful  to  current  ongoing  projects  when  it  is  treated  with  more 
modem  methods  of  data  analysis.  In  specific  the  report  does  the  following 
things: 


1.  It  transforms  old  data  suspected  or  shown  to  have  relevance  to  current 
projects  into  graphical  func*.lons  which  represent  the  trend  of  this  data. 

2.  It  shows  that  these  graphical  functions  can  be  used  tn  prcv?de  new 
information  capable  of  being  used  in  further  a-iaXysas. 

3.  This  new  information  is  shown  to  form  new  tables  which  can  be  analyzed 
with  methods  derived  from  response  surface  methodology  to  produce  a  more 
generalized  model  capable  of  more  broadly  representing  all  of  the  data  trends. 

4.  This  model  nay  be  used  to  make  predictions  which  would  be  difficult  to 
do  rlth  Che  original  data. 

5*  The  model  -  Chot^h  eiqiirical  in  nature  -  may  be  extremely  helpful  in 
reducing  the  level  of  experimentation  required  to  develop  the  necessary  data 
bases  needed  in  ongoing  casks  or  projects. 

6.  The  model  may  also  suggest  more  neaningfui  experiments  and  or  reduced 
experimentacion  to  cover  the  goals  and  objectives  of  current  ongoing  casks. 

Such  reductions  in  experimentation  translate  into  reduced  costs  and  earlier 
accomplishment;  of  the 

The  report  is  an  outline  of  what  can  be  done  with  old  data  buried  in  the 
licecacuret  or  even  in  old  notebooks  where  there  is  a  degree  of  reasonable 
completeness.  Under  the  Data  Analysis  and  Methodology  Section  it  shows  chat 
Che  old  data  can  be  transforxMd  into  sim^och  functions.  In  the  Results  and 
Conclusions  Secclon  it  shows  that  these  smooch  functions  can  be  incorporated 
into  a  generalized  model  which  can  model  all  Che  seemingly  diverse  original 
data.  It  is  also  shown  in  this  section  chat  three  dimensional  graphs  generated 
from  the  function  representing  the  model  reveal  what  the  model  looks  like  and 
give  a  strong  feeling  for  the  influence  of  the  impoitant  parameters  on  the 
principal  property  represented  by  the  model. 

With  respect  to  the  specific  illustrative  example  of  this  report  for  the 
interaction  of  mustard  gas  with  paint  the  following  results  and  conclusions 
were  obtained: 

a)  Models  were  developed  which  represented  mustard  paint  interactions 
over  a  range  of  temperatures  and  wind  speeds. 


b)  These  models  represenced  the  effect  of  temperature  and  wind  speed  on 
the  mustard  evaporation  time  and  the  sorption  power  of  the  paint  for  mustard 
liquid. 

c)  The  sorption  power  model  shows  that  fresh  paint  absorbs  mustard  more 
readily  than  aged  paint. 

d)  The  sorption  power  model  also  shows  chat  temperature  promotes 
absorption  of  mustard  into  paint  while  wind  speed  lessens  absorption  into  the 
paint. 

e)  The  evaporation  models  show  chat  both  temperature  and  wind  speed 
promote  effective  evaporation  of  mustard  and  that  these  evaporation  times  are 
reasonably  cracked  by  the  model  over  the  experimental  range  of  temperature  and 
wind  speed. 
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THE  MUSTARD  GAS-PAINT  SYSTEM  REVISITED 


INTRODUCTION 

CW  decontamlnatioa  is  a  costly  and  time-consuming  operation  which  Is  both 
manpower-intensive  and  manpower-constrictive.  The  need  for  decontaalnation  may 
be  properly  addressed  when  the  type  and  concentration  of  the  deposiced  chemical 
agent  are  known.  Weathering  becomes  a  suitable  alternative  to  decontamination 
when  the  environmental  parameters  of  temperature.  Incident  radiant  energy  and 
wind  speed  are  shovm  to  promote  favorcble  evaporation  of  the  deposiced  chemical 
agent  with  respect  to  operational  considerations. 

Mustard  gas  (HD)  is  a  favorite  chemical  agent  for  use  as  a  laboratory 
model  of  a  chemical  agent.  While  it  lacks  the  killing  power  of  Che  exotic 
nerve  agents  by  at  least  an  order  of  magnitude,  Ics  persistence  and  abllicy  to 
cause  disfiguring  life-long  injuries  and  incapacitation  to  survivors  is 
well-known.  It  also  has  the  advantage  of  tying  up  signlficanc  numbers  of 
personnel  to  administer  to  Che  injured  individuals  as  opposed  to  merely  setting 
aside  Che  dead  for  burial.  It  (HD)  also  is  one  of  the  more  intractable  agents 
to  decontaminate  in  addition  to  being  relatively  persistent.  Mustard,  in  spite 
of  ics  long  history,  remains  the  agent  of  choice  in  practical  studies  of 
deccntamination  and  weathering  effects. 

One  of  the  early  laboratory  weathering  studies  on  mustard  gas  was  done  at 
this  Laboratory  in  1961  (1)  by  A.  Scamulis.  This  report  looked  at  the  evapora¬ 
tion  of  BO  from  both  fresh  and  aged  painted  Navy  marine  surfaces  as  a  function 
of  temperature  and  wind  speed  under  controlled  laboratory  conditions.  In  this 
report,  Scamulis  treated  Che  evaporation  process  as  being  divided  into  two 
characteristic  regions:  (a)  the  exposed  liquid  surface  region  and  (b)  the 
sorbed  liquid  (into  the  paint)  region.  Most  importantly,  Stamulis  provided 
data  cables  for  his  controlled  experiments  which  are  most  worthy  of  being  re¬ 
examined  in  Che  light  of  new  methcdoiogies. 

This  report  will  revisit  this  data  and  accempt  to  develop  relations  which 
may  be  useful  for  chose  concerned  with  the  weathering  properties  of  the 
mustard-paint  system.  If  the  results  presented  In  this  report  have  some  aerie 
or  application  or  if  the  methods  used  will  be  useful  in  making  practical 
guesses  or  judgments  with  respect  to  weathering  decisions,  then  this  report 
can  be  justified.  An  important  objective  of  this  report  is  to  be  able  to 
predict  the  percent  of  mustard  evaporated  from  paint  surfaces  as  a  function  of 
temperature,  wind  speed  and  time  from  a  minimum  of  experimental  data.  In  order 
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CO  achieve  this  objective  a  conbinaCion  of  analytic  methods  will  be  used  which 
will  (I)  reduce  Scaaulis's  evaporation  curve  results  to  a  series  of  two  para¬ 
meter  equations  and  (2)  develop  an  empirical  equation  which  will  show  the 
effect  of  Ce]q)erature  and  wind  speed  using  the  evaporation  lialf  life  parameter 
developed  from  the  two  parameter  evaporation  curve  plots  (Michaelson-Menten 
plots) . 

This  report  also  features  the  use  of  response  surface  methodology  (RSM) 
methods  to  extract  relationships  from  the  experimental  data  array.  In  order  to 
make  progress  rapidly  in  the  use  of  a  new  method  or  technique  It  is  necessary 
to  (1)  gain  a  knowledge  of  how  to  apply  the  given  method  or  technique  to  the 
problem  of  interest  and  (2)  gain  an  understanding  of  the  principles  arid  theory 
on  which  the  method/technique  is  based.  It  is  the  writer's  contention  chat  the 
two  are  to  some  degree  auctially  exclusive  in  that  a  thorough  understanding  of 
the  principles  on  which  a  particular  technique  is  based  is  not  absolutely 
necessary  to  its  initial  application.  This  report  will  not  attempt  to  explain 
the  statistical  theory  behind  the  use  of  RSM  but»  rather >  quickly  expose  the 
reader  to  the  tools  by  which  he  may  organlze/plan/interpret  his  experiments, 
rhe  reference  section  contains  information  on  the  statistical  theory  behind  the 
development  of  R$K  ^2,3).  Upon  completion  of  these  operations*  the  reader  will 
be  left  with  a  reasonable  view  of  mustard  evaporation  as  a  function  of 
environmental  parameters  and  paint  type. 


DATA  A2^TSIS  AND  METHODOLOGY 
Mlchfi^ls'ir-41enten  Analysis 

i*relimiiiary  analysis  of  Staatulls's  data  showed  chat  the  sorbed  agent  dora 
could  be  reasonably  well  represented  by  a  hyperbolic  equation  of  the  type 

2EV  -  (c/(t  +  K))*  2EV(maX) 

where  K  and  ZEV(MAX)  are  constants  Co  be  determined  by  the  data  fit*  t  •  time 
in  hours  and  ZEV  *  percent  mustard  evaporated  at  time  t.  The  constant  K 
represents  Che  half  life  of  the  evapora’’tco  process*  that  is  the  point  in  time 
at  which  one  half  of  the  deposited  c«.<naical  agent  has  disappeared.  The 
constant  ZEV(MAX}  is  ch«  pe/cent  of  agent  evaporated  av  infinite  time  which 
ideally  should  be  IGOZ  in  vclue. 

There  are  three  ways  ta  plot  the  experimental  data  ro  such  an  equation  in 
order  to  obtain  these  parameters  (4).  A  simple  laast  squares  fit  computer 
program  adapted  from  a  Xfteriture  Fortran  program  been  used  to  show  the 
difference  between  chess  plot  results.  The  .jstt.cds  are  (1)  the 

Lineweavec-Burk  method,  (2)  Che  HichaeXsoa'’dence'.t  method  and  (3)  the 
Eadie-Hofscee  method.  The  Lineweaver-Burk  nechoo  plots  i/Z£V  vs  i/c.  The 
slope  is  ZR/Z£V(MAX}  while  the  intercept  with  Che  y  axis  is  ;/ZEV(MAX)  and  with 
th(>  X  aria  is  -l/R.  A  problem  with  this  XKSthod  is  that  data  points  with  small 
values  which  are  inaccurate  are  too  strcngly  weighted  when  reciprocals  are 
plotted.  The  Michaelson-Menten  method  plots  c/Z£V  vs  c.  The  slope  of  this 
plot  is  1/Z£V(MAX)  and  the  intercept  with  Che  y  axis  is  R/ZEV(MAK).  In  the 
Eadie-Hofscee  method  ZEV  is  plotted  against  ZEV/c.  The  slope  is  -K  and  the 
intercept  with  the  y  axis  is  ZEV (MAX)  and  with  the  x  axis  is  ZEV(ma.S)/K.  The 
method  selected  In  this  report  was  the  Michaelis-Menteo  method  since  present 
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experience  shows  it  to  give  the  most  consistent  and  reliable  results  over  the 
observed  range  of  data  points.  For  example.  Table  I  shows  the  results  obtained 
by  least  squares  using  the  three  methods  for  two  quite  different  data  files. 

The  Lineveaver-Burke  Method  invariably  gives  the  lowest  values  of  variance  and 
standard  deviation  but  the  values  of  the  parameters  are  sometimes  far  afield  of 
the  expected  values.  The  Eadie-Hofstee  Method  gives  values  which  are  more  in 
line  with  expectations  but  the  variance  and  standard  deviation  are  almost 
always  significantly  greater  chan  for  the  Michaelson-Menten  Method.  Figure  1 
is  an  example  of  the  data  fit  to  the  Michaelis'Mencen  relation  where  K  was 
found  to  be  0.3616  and  the  value  of  Z£V(HAX)  was  found  to  be  lOOZ.  Plots  of 
Stamulis's  data  at  various  temperatures  and  windspeeds  were  done  by  computer 
and  reasonable  conformity  to  the  relation  was  found.  The  values  of  K  were  sec 
aside  for  each  condition  of  temperature  and  wind  speed. 

Response  Surface  Methodology 

A.  useful  RSM  design  adaptable  to  Stamulis's  data  is  the  face  or 
side->cencered  design  (we  will  not  attempt  the  pros  and  cons  of  one  design  vs 
another  in  this  report) .  This  design  produces  a  polynomial  which  contains  a 
constant  term,  linear  terms,  quadratic  terms  and  cross-product  terms  and  is  of 
the  form: 


^  -  Bo  +  i;  BiXi  +  £  +  f  S  Bi^XiXj 

i»l  i*l  i*l  j>i 


In  order  to  develop  the  coefficients  of  the  above  polynomial  it  is  necessary  to 
run  experiments  in  the  independent  variables  (x’s)  and  observe  Che  dependent 
variables  (y's).  To  successfully  accomplish  this  goal  requires  an  experimental 
'recipe*  which,  when  used  with  an  appropriate  array,  will  allow  computer 
solution  of  the  desired  coefficients.  The  coefficients  obtained  are  a  least 
squares-  like  fit  of  the  data  to  the  polynomial.  It  should  again  be  emphasized 
that  the  polynomial  obtained  is  not  a  representation  of  the  'natural  law'  which 
is  at  work.  It  will,  however,  represent  the  variation  of  the  response 
variables  (y's)  to  the  independent  variables  (x's)  perhaps  as  veil  as  the 
natural  lav  Itself  over  the  same  range  of  independent  v^iriables  from  which  the 
polynomial  was  obtained  were  the  natural  law  adequately  known. 

In  the  case  of  2  variables,  the  number  of  experiments  for  single  runs 
without  duplication  should  be  9  (4  for  each  comer  of  the  square  plus  4  center 
points  in  each  side  plus  the  center  of  the  square).  The  data  array  required 
for  a  polynomial  containing  Z  independent  variables  with  interactions  is  shown 
in  Table  2.  Note  that  the  x*s  are  given  in  coded  form.  This  is  an  aid  in  both 
conforming  to  the  experimental  design  ard  in  the  rational  solution  of  the 
equations  set  up  for  determining  the  coef ficZ.ents.  The  2  variable  situation  is 
probably  one  of  the  most  common  encountered  in  the  laboratory  and  its  inclusion 
here  will  make  a  reasonable  example.  The  half  matrix  representing  the  normal 
equations  is  shown  in  Table  3.  The  specific  relations  for  Bl,  B2  and  B12  are 
shown  together  with  the  full  reduced  matrix  which  must  be  inverted  to  solve  for 
BO,  BU  and  B22. 
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One  o£  the  principal  thrusts  behind  the  intent  to  produce  a  statistical 
equation  is  the  development  of  a  function  by  means  of  which  surface  plots 
and/or  contour  maps  may  be  generated.  Such  maps  will  greatly  aid  the 
experimenter  in  obtaining  a  grand  view  of  his  experimental  results  over  the 
range  of  interest  and  alsot  perhaps,  to  make  some  limited  predictions.  Figure 
2,  then,  stnum^rizes  the  operations  and  goals  described  in  this  section. 
Examples  of  RSM  methods  and  software  used  by  the  author  relevant  to  this  paper 
may  be  found  In  earlier  reports  (5, 6, 7, 8). 


III.  RESULTS 

Tabulation  of  Data  Hatrix  (RSM  Analysis) : 

Table  4  is  a  summary  of  Stamulls's  evaporation  data  taken  at  a  variety  of 
conditions  and  on  several  types  of  Havy  paints  both  fresh  and  aged.  The 
conditions  ranged  from  30*C  to  S0*C  in  temperature  and  from  0.1  to  16.1  miles 
per  hour  in  wind  speed.  The  paints  used  were  of  the  alkyd  resin  type,  both 
fresh  and  aged  (eight  years  old).  Also  included  in  the  Table  are  data  computed 
from  Che  Michaelson-Menten  plots  which  were  recorded  as  corresponding  values 
of  R  and  2EV(MAZ)  for  each  condition  of  temperature  and  wind  velocity.  Table  4 
provides  the  basis  for  the  development  of  Tables  3  and  6  which  contain  selected 
and  interpolated  values  of  the  independent  and  dependent  (response)  variables 
to  conform  to  the  orthogonal  side  centered  statistical  design  mentioned 
earlier.  Table  5  summarizes  the  extracted  data  for  the  aged  paint  experiments 
while  Table  6  presents  Che  mustard  evaporation  data  from  freshly  painted 
surfaces.  It  should  be  noted  that  Staaulis's  data  mostly  conformed  to  the  side 
centered  statistical  design  with  only  minimal  interpolations  being  needed  to 
arrive  at  the  $.1  mph  condition  required  for  the  orthogonality  condition 
instead  of  Che  6.1  xoph  data  documented  in  Scamulis's  report.  In  the  design  it 
should  be  noted  that  the  dependent  variable  is  taken  as  the  logarithm  of  the 
half  life  of  Che  given  muscard-^painc  system.  Such  a  dependent  variable 
transformation  allows  a  much  closer  correlation  with  Che  independent  variables 
of  the  experimental  design.  Also  included  in  Tables  5  and  6  are  data  on 
Scamulis's  "sorbing  power  coefficient",  Ksp.  These  data  were  also  entered  for 
a  total  of  4  separate  equation  deteraxnaclons.  Staaulis  defines  his  sorption 
power  coefficient  as  "the  ease  with  which  a  given  paint  film  absorbs  mustard 
gas  under  given  evaporation  conditions".  Ko  cransfomation  of  the  response 
variable  was  required  for  the  sorption  power  coefficients. 

Development  of  Descriptive  Equations: 

Table  7  summarizes  the  coefficients  obtained  by  solution  of  the  four 
systems  of  normal  equations  represented  by  the  matrices  of  Table  2  and  3 
together  with  the  response  variable  data  of  Tables  5  and  6.  These  coefficients 
are  substituted  as  appropriate  in  the  given  descriptive  equation.  As  mentioned 
in  the  earlier  section,  these  equations  may  now  be  used  to  (1)  generate  three 
dimensional  surfaces  and/or  contour  maps  and  (2)  make  predictions  of  experi¬ 
mental  behavior  covered  by  the  experimental  design.  It  must  be  emphasized  chat 
experimental  designs  of  this  type  cannet  be  used  to  make  extrapolations  beyond 
the  region  enclosed  by  the  design. 


Predictive  Results  vs  Actual  Data: 


Table  8  presents  the  results  of  computing  the  percent  mustard  evaporoted 
using  K  derived  from  the  descriptive  two  variable  equations  of  the  previous 
section*  one  equation  representing  aged  paint  and  the  ether  the  fresh  paint 
relation.  Only  the  R  value  was  used  in  the  computations  of  the  percent 
evaporated  as  ZEV(MAX)  was  set  to  lOOZ.  Thus*  coded  values  of  temperature  and 
vindspeed  were  substituted  Into  the  descriptive  equation  to  generate  the  K 
value.  The  K  value  was  substituted  Into  the  Michaelson-Henten  hyperbolic 
relation  assuming  ZEV(MAX)  *  lOOZ  and  the  data  values  calculated  for  the  time 
intervals  Indicated  In  Table  8.  In  spite  of  the  fact  that  some  of  Stamulis's 
data  was  rather  scanty  and  frankly  incomplete*  the  agreement  between  the 
predicted  results  and  the  original  data  is  reasonable  in  most  cases  and  gives 
some  confidence  in  looking  at  the  general  appearance  of  the  three  dimensional 
surface  function  generated  by  the  descriptive  equation. 

Evaporation  and  the  Paint  Sorption  Power  Coefficient: 

Stamulls  has  defined  a  paint  sorption  power  coefficient  based  on  a 
comparison  of  his  experimental  evaporation  data  with  his  "dry  line"  curve.  No 
attempt  will  be  made  here  to  go  into  the  details  of  his  "dry  line"  concept  but 
rather  to  point  out  that  a  straight  forward  analysis  of  the  evaporation  from 
paint  films  as  opposed  to  evaporation  from  bare  metal  surfaces  leads  more 
clearly  to  the  same  results  for  the  value  of  paint  sorption  power  coefficients 
as  he  obtained  in  his  report  with  no  graphical  geometry  required  in  the 
analysis. 

figure  3  is  a  diagram  of  chemical  agent  evaporation  from  a  bare  metal  as 
opposed  to  a  paint  film.  The  time  required  to  evaporate  lOOZ  of  the  agent  from 
the  bare  metal*  t'*  is  used  to  estimate  tha  value  of  ZEV* .  It  is  clear  that 

ZEV*  ■  (t'/(t*  +  K))100  or  ZRV/lOO  -  tV(t*  +  K) 

If  the  sorption  coefficient  is  defined  as  Ksp  *  (100  •  Z£V')/100*  Ksp  may  be 
reduced  to  Ksp  1  -  ZEV'/IOO.  Upon  substitution  into  the  above  equation  one 
obtains 

Ksp  -  K/(t'  +  K) 

Comparison  of  sorption  power  coefficients  calculated  directly  from 
mustard-paint  film  half  lives  showed  an  agreement  of  +-  lOZ  with  Stamulis's 
original  Ksp  values  without  the  use  of  the  any  graphical  geometry  for  the 
computation. 

Three  Dimensional  Graphics: 

Figures  4  through  8  show  the  power  of  three  dimensional  graphics  giving  a 
broad  visual  view  of  exactly  how  a  function  is  behaving  over  the  desired  range 
of  interest.  Figures  4  and  S  compare  the  behavior  of  mustard  evaporating  from 
aged  paint  with  evaporation  from  fresh  paint.  The  x-axis  represents  the  effect 
of  temperature  and  the  y-axis  the  effect  of  wind  speed  while  the  z-a\is  is  the 
normalized  effect  of  these  two  independent  variables  on  logarithm  of  the 
mustard  half-life.  Figure  6  is  a  plot  of  the  difference  between  the  logarithms 
of  the  half-lives  of  the  aged  anc  fresh  painted  surfaces  vs  temperature  and 


wind  speed.  Vhere  values  on  Che  z-axis  are  low  ic  indicates  chat  curve 
differences  are  small  and  vice  versa.  It  would  appear  chac  the  paints  (aged  vs 
fresh)  show  the  greatest  differences  at  higher  wind  speeds  since  the  difference 
in  the  logarithms  of  the  half-lives  is  highest  at  the  greater  wind  speeds. 

Figures  7  and  8  compare  the  paint  sorption  power  coefficients  as  affected 
by  temperature  and  wind  speed  for  aged  paint  vs  fresh  paint.  Inspection  of 
Table  7  and  comparison  of  the  figures  show  vividly  how  the  fresh  paint  more 
effectively  absorbs  the  mustard  as  against  the  aged  paint  and  how  an  increase 
in  temperature  promotes  increased  sorption  of  the  muscara  gas  by  the  paint  in 
both  cases.  As  expected.  Increased  wind  speed  decreases  the  sorption  power 
coefficient  for  both  aged  and  fresh  paints.  In  the  case  of  these  alkyd  type 
paints  while  temperature  increases  the  race  of  evaporation,  ic  also  promotes 
the  penetration  of  the  chemical  agent  into  the  paint.  Moreover,  the  aging  of 
alkyd  type  paints  appears  to  Inhibit  the  penetration  of  the  cnemlcal  agent  into 
the  paint  film. 


IV.  CONCLUSIONS 

a)  Models  can  be  developed  which  represent  the  dependence  of  important 
system  properties  over  a  range  of  independent  variables  when  properly  selected 
archival  data  is  used  as  a  basis  for  these  models. 

b)  The  exemplary  models  of  this  report  represent  the  effect  of  temperature 
and  wind  speed  on  the  mustard  evaporation  time  and  the  sorption  power  of  the 
paint  for  mustard  liquid. 

c)  The  sorption  power  nodel  shows  chac  fresh  paint  absorbs  mustard  more 
readily  chan  aged  paint. 

d)  The  sorption  power  model  also  shows  that  temperature  promotes  absorption 
of  mustard  into  paint  while  wind  spt^d  lessens  absorption  into  the  paint. 

e)  The  evaporation  models  show  chat  both  temperature  and  wind  speed  promote 
effective  evaporation  of  mustard  and  that  these  evaporation  times  are 
reasonably  veil  cracked  by  the  model  over  the  experimental  range  of  temperature 
and  wind  speed. 


V.  RECOMMENDATIONS 

A  large  database  of  investigated  chemical  agent  systems  lies  in  the  variety 
of  government  reports  published  over  the  years.  While  these  reports  are 
accessible  through  DOD  and  NTIS  services  etc.  many  reports  such  as  the  one 
examined  here  may  have  new  value  if  re-examined  within  the  context  of  new 
analytic  metheds  or  methodologies.  Since  the  costs  of  live  agent  tests  and 
research  in  general  have  soared  in  recent  years,  it  might  be  useful  to 
critically  examine  the  older  reports  in  a  systematic  manner  using  the  newer 
analytic  methods  presented  here  to  reduce  some  of  the  expensive  and  time 
consuming  live  agent  testing  to  which  the  armed  services  are  committed.  In 
specific  and  to  the  point: 
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1.  That  a  ta^k  be  set  up  to  critically  review  the  content  of  WVll  and  post 
Wtfll  chemical  agent  work  to  determine  the  merit  of  this  work  against  the 
requirements  and  planned  research  addressed  by  current  ongoing  or  planned  CBR 
tasks. 


2.  That  the  more  useful  and  relevant  of  these  reports  be  reformatted 
using  methods  similar  to  the  ones  used  in  this  report  in  order  to  both  extract 
the  maximum  amount  of  information  and  reduce  the  cost  of  additional  live  agent 
testing-  especially  if  it  is  established  that  similar  work  has  already  been 
accomplished  in  this  earlier  literature  relevant  to  these  CBR  casks. 
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Table  1:  Lease  Squares  Data  Plots  for  the  Three  Methods 
HD.IO.STAM.ETK  HD.33.STAM.ETK 


s 

V 

s 

V 

1.5 

30.1 

1 

11 

2 

35.6 

2 

22.5 

3 

48 

3 

35.5 

4 

59.5 

5 

62.4 

5 

68.6 

6 

74.9 

6 

73.1 

7 

80 

24 

90.5 

24 

87 

48 

94.2 

48 

88.2 

LINEWEAVER-BURK  METHOD 

1 

as 

M 

-BURK  METHOD 

ZEV(MAX)  - 

116.364527 

ZEV(MAX)  » 

333.171971 

K  -  4.27860016 

K  »  28.2869931 

VAP'.-NCE  - 

2.2535524E-09 

VARIANCE  • 

1.22701418E-03 

STND  DEV  - 

4.74715957E-05 

STN  DEV  -  i 

1.22701418E-03 

MICHEALIS-MEMTEH  METHOD 

ZEV(MAX)  -  100.493319 
K  «  2.95253882 
VARIAHCE  •  5.09142197E-06 
STS  DEV  «  2.25641795E-03 


MICHEALIS-MESTES  METHOD 

ZEVfMAX)  •  98.6486888 
K  -  4.42648385 
VARIANCE  -  2.14812686E-04 
STSD  DEV  •  .0146564895 


EADIE-HOFSTEE  .METHOD 

ZEV(MAX)  •  106.271227 
K  -  3.49248494 
VARIANCE  -  6.17870316 
STN  DEV  -  2.48569973 


EADIE-HOFSTEE  METHOD 

ZEV(MAX)  •  94.0696557 
K  -  3.3328406 
VARIANCE  »  278.916323 
STS  DEV  -  16.7007881 
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Table  2.  Side-Centered  Two  Variable  Design 


Design  Matrix  Graphic  Representation 


H 

Xl 

X2 

Y 

j 

1 

1 

1 

Itl 

X 

X 

X 

2 

i 

-1 

*2 

Comers  j 

3 

-1 

1 

'13 

\ 

u 

-1 

-1 

V4 

X2 

X 

X 

X 

5 

1 

0 

ifs 

Side  of  Face 

6 

-1 

0 

*6 

Centered  Pts  1 

7 

0 

1 

If? 

X 

X 

X 

8 

0 

-1 

TfS 

Center  Point 

9 

0 

0 

lf9 

Xl 

Table  3.  Solution  of  Kormal  Equations 

u 

EXl 

ZX2 

EXi^ 

2x2^ 

1X1X2 

ZY 

rxiZ 

1X1X2 

ZXi3 

1X1X2^ 

EXi2X2 

ZXiY 

EXi2X2 

EX23 

EX1X22 

i:X2Y 

EXi* 

rxi2x22 

2Xi3x2 

:xpY 

2X2'* 

ZX1X23 

r:<23Y 

IXi^X,^ 

i:XlX2Y 

9 

0 

0 

6 

6 

0 

r; 

6 

0 

0 

0 

0 

:xiY 

6 

0 

0 

0 

LX2Y 

4p 

6 

4 

0 

IXi-Y 

6 

0 

ZXi^Y 

4 

31X2Y 

Bi  - 

EXiY/6 

Note: 

B2  - 

nC2Y/6 

Relations  for 

Linear  and  Intereaction 

Coefficients 

®12  ■ 

■  rXiX2Y/4 

9  6  6  n  Note: 

664  Invert  Matrix  to  Solve  for  and  Second  Order 

6  4  6  IX.2^^  Coeficients 
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Table  5.  Eight  Year  Paint  Data  Used  in  Design 


STAHULIS 


EXP. 

NO.  TEMP 

*1 

WIND 

*2 

Y 

tl/2 

!.ogtl/2 

kjp 

18 

50 

1 

16.1 

1 

Y(l) 

.0451 

-1.346 

.23 

15 

50 

1 

0.1 

-1 

Y(2) 

.385 

-.4145 

.23 

10 

50 

1 

16.1 

1 

Y(3) 

.0521 

-1.283 

.08 

07 

30 

-1 

0.1 

-1 

Y(4) 

1.640 

.2148 

.15 

* 

50 

1 

8.1 

0 

Y(5) 

(.0720) 

-1.143 

(.18)est. 

* 

30 

-1 

8.1 

0 

Y(6) 

(.0808) 

-1.092 

(.12)esc. 

14 

40 

0 

16.1 

1 

Y(7) 

.0505 

-1.296 

.13 

11 

40 

0 

0.1 

-1 

Y(8) 

1.807 

.2570 

.35 

*  40  0  8.1  0 

*  •  interpolated  row  values  from  data 

Y(9) 

(.0408) 

-1.389 

(.13)est. 

STAMOLIS 
EXP.  NO. 

TEMP 

Xl 

Table  6.  Fresh  Paint  Data 

WIND  X2  Y 

Used  in  Design 

Cl/2  Logtl/2 

18 

50 

1 

16.1 

1 

Y(l) 

.1046 

-.9805 

.23 

27 

50 

1 

0.1 

-1 

Y(2) 

.3650 

-.0530 

.53 

22 

30 

-1 

16.1 

1 

Y(3) 

.3394 

-.4693 

.39 

19 

30 

-1 

0.1 

-1 

Y(4) 

2.952 

.4701 

.33 

* 

50 

1 

8.1 

0 

Y(5) 

(.2424) 

-.6155 

(.48)est. 

* 

30 

-1 

8.2 

0 

Y(5) 

.2994 

-.5237 

(.23)esc. 

26 

40 

0 

16.1 

1 

Y(7) 

,2000 

-.6990 

.32 

40 

0 

0 

O.l 

-1 

Y(8) 

.8658 

-.0626 

.80 

* 

40 

0 

8.1 

0 

Y(9) 

(.1258) 

-.9003 

(.40)esC. 

*  >  interpolated 

row 

values  froa  data 
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Table  7.  SuoDary  of  Equation  Coefficients 


COEFFICIENT 

AGED  PAINT 

Y-LogK  Y  -  Ksp 

tRESH 
Y  «  LogK 

PAINT 

V 

^sp 

• 

Bo 

-1.2579 

0.1550 

-0.8492 

0.4343 

- 

Bl 

-0.1239 

0.04833 

-0.1684 

0.0650 

B2 

-0.6637 

-0.04833 

-0.4366 

-0.i367 

Bll 

0.00914 

0.03000 

0.2284 

-0.1136 

*22 

0.6072 

0.0600 

0.4172 

0.09U 

*12 

0.1416 

0.0175 

-0.0260 

-0.1150 

MODEL  EQUATION:  Y  -  Bq  +  BiXi  +  83X2  +  BuXi^  +  832X2^  +  812^1X2 


X^  *  Coded  Tesperature 

X2  *  Coded  Wind  Speed 

xi  -  f:!* 

X,  -  ws  -  8.1 

10 

8 

12 


Trtbltt  8.  Actual  Data  va.  Prcdlctlva  Naaulta  (Jn  Parinthaaua) 
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FIGURE  I.  Evnporatlon  of  Miiii(ard  from  Fresh  Faint:  Comparison  of 
Equation  and  Experimental  Data  Points. 


FURTHER  DATA  ANALYSIS  j 
1.  FACTOR  INTERACTIONS  i 
I  2.  PREDICTIONS 
I  3.  OPTIMIZATION 


Figure  2.  Treatoenc  of  Experimencal  Daca  and  Analysis 


Klguic*  J.  r.viiporiition  of  Miiaiord  From  NcUl  vs.  AlUyd  Kusln  Faint: 
DutuiminaClon  ol  t*  aiiJ  XJ’.V’. 


EVAPORATION  HALF  LIFE 
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Figtiru  4.  Kvaporntlon  ^fodel  ok  Miihtard  From  Agud  Faint  an  a  Function 
of  Tumpuraturu  and  Wind  Spuud. 
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EVAPORATION 
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Figure  6.  Graplilcul  Siibtruccion  of  Aged  Pulut 


SORPTION  POWER  COEFFICIENT 


Function  of  Tumperntiire  niul  Wind  Speed  for  Aged  Paint 


POWER  COEFFICIENT 


Mgtiru  8.  Sorpcloi)  Po\/er  (\)effJciu»L  a  riiiiclion  of  TuinpuraCuru  am)  Wind  Speed  for  Ftcblt  t’«iint  Mixlel. 


